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ABSTRACT 



Context. The knowledge of the present-day mass function of young clusters and the mass of their coolest substellar members is es- 
sential to clarify the brown dwarf formation mechanism, which still remains a matter of debate. 
Aims. We searched for isolated planetary-mass T-dwarfs in the ~3 Myr old Serpens Core cluster. 

Methods. We performed a deep imaging survey of the central part of this cluster using the WIRCam camera at the CFHT. Observations 
were performed through the narrow-band CH 4 off and CH 4 on filters, to identify young T-dwarfs from their 1.6yum methane absorption 
bands, and the broad-band JHKs filters, to better characterize the selected candidates. We complemented our WIRCam photometry 
with optical imaging data from MegaCam at CFHT and Suprime-Cam at the Subaru telescope and mid-infrared flux measurements 
from the Spitzer "core to disk" (c2d) Legacy Survey. 

Results. We report four faint T-dwarf candidates in the direction of the Serpens Core with CH 4 on-CH 4 off above 0.2 mag, esti- 
mated visual extinction in the range 1-9 mag and spectral type in the range T1-T5 based on their dereddened CH 4 on-CH 4 off colors. 
Comparisons with T-dwarf spectral models and optical to mid-infrared color-color and color-magnitude diagrams, indicate that two of 
our candidates (ID 1 and 2) are background contaminants (most likely heavily reddened low-redshift quasars). The properties of the 
other two candidates (ID 3 and 4) are consistent with them being young members of the Serpens Core cluster, although our analysis 
can not be considered conclusive. In particular, ID 3 may also be a foreground T-dwarf. It is detected by the Spitzer c2d survey but 
only flux upper limits are available above 5.8 /jm and, hence, we can not assess the presence of a possible disk around this object. 
However, it presents some similarities with other young T-dwarf candidates (S Ori 70 in the a Orionis cluster and CFHT J0344+3206 
in the direction of IC 348). If ID 3 and 4 belong to Serpens, they would have a mass of a few Jupiter masses and would be amongst 
the youngest, lowest mass objects detected in a star-forming region so far. 
Conclusions. 

Key words, stars: formation - stars: low-mass, stars: brown dwarfs, ISM: clouds, ISM: individual objects: Serpens Core 



1. Introduction 

Understanding the brown dwarf (BD) formation mechanism is a 
key point to assess what determines and which is the minimum 
mass for star formation. Although two decades h ave passed since 
the fi r st unambiguo us observations of BDs (Nakaiima et al. 
Il995t iRebolo et all 1 19951) . this issue is still under debate. 
According to the star formation theory, the minimum mass for 
star formation is set by the so called mass fragmentation limit, 
i.e. the mass at which one object cannot contract further because 
the radiated heat preven ts it to collapse further. This limit is of 
the order of a 5-10My dLow & Lvnden-Belllll976t lReeslll976l) 
and may be lower (~1 My) in the pr esence of magnetic fields 
dBossll200lUBovd & Whitworthll2005T) . Observationally, several 
young star forming regions have been probed to faint magnitudes 
at both optical and near-infrared wavelengths to look for a pos- 
sible cut-off at the low-mass end of the mass function and the re- 
sults significantly vary from region to region. A possible cut-off 
below ~6 My up to 20 My have been reported in NGC 1333 (e.g., 
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IScholz ef al. 2009, 2012). Moreover, there is no universal agree- 
ment on the behavior of the mas s function over this c ut-off. A 
search for T-dwarfs in IC 348 by Burgess et al. ( 2009) suggests 
that a n extrapolation of the log-normal mass function jChabrierf 
2003) may hold in the substellar and planetary-mass regimes. A 
similar study in the Upper Sco associ ation favors a tur n-down 
in the m ass function below 10 My dLodieu et alj|201 it) . In the 
Pleiades Casewell et al. (2007) reported several L/T-type candi- 
dates with masses as low as 10 My, implying that the slope of 
the mass function in this mass regime agrees, within the uncer- 
tainties, with the values inferred from earlier studies at higher 
masses (ICasewell et al.ll2007t bobbie et al.l 12002b lLodieu et al.l 
120071: iMoraux et alj2003l) . Finally. iKirkpatrick et al.ld2012l) find 
tantalizing hints that the number of BDs in the field continues to 
rise from late-T to early-Y. 

To address these questions, a large observing key-program 
(P.I. J. Bouvier) was conducted at the Canada France Hawaii 
Telescope (CFHT) aimed at characterizing the stellar population 
of nearby star forming regions using the Wide Field IR Camera 
(WIRCam). One of the main objectives of this program is to 
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look for the lowest mass members of these clusters and inves- 
tigate the cut-off at the low mass end of the mass function, us- 
ing a technique bas ed on narrow-band methane imaging (e.g., 
iBurgess et al]l2009t) . 

In this paper, we report the results of this observing pro- 
gram in the Serpens Core cluster, specifically focusing on the 
methane filters imaging to search for young T-dwarfs can- 
didates. The Serpens Core region is an example of a very 
young, deeply embedded (Ay =s40 mag) cluster , containing 
a hig h percentage of protostars dDavis et al.l 1 19991: lEiroa et alj 
2006: lHarvev et alj l2007t iKaas et all l2004t iTesti et alj 12000: 
Winston et al] |2007h . With an age betw een 2 and ~6 Myr 
(Oliveira et al. 20091: IWinston et al] 20091) and a distance esti- 
mated between 2 60+37 pc dStraizvs et al]|2003l) and 415±5 pc 
(Dzi b et aDl2010h . the Serpens cloud core is one of the nearest 
regions of clustered star formation to the Sun. Therefore, it is 
an excellent candidate for study as it is close enough to both re- 
solve the individual members and detect the lowest mass mem- 
bers to below the hydrogen-burning limit. Through out this paper, 
we assume a typical age of 3 Myr (IWinston et al. 2009) for the 
Serpens core cluste r and the minimum distance of 260+37 pc 
dStraizvs etalJl200l R The region is also particularly interest- 
ing for BD studies because only a few dedicated searches have 
been reported in t he literature. The fi rst young BD in Serpens 
was discovered bv lLodieu et alj d2002l) and the cluster substellar 
population count so far about 3 spectrosc opically confirmed BDs 
(Lodie u et al ] l2002HShirono et al.ll2011l) and 40 BD candidates, 
a considerable fraction of which still surrounded by prominent 
accretion disks (lEiroa et alJl2006h . 

Canonical approaches on BP studies and related issues (see, 
e.g., iJavawardhana et all 120031) rely on the identification and 
spectral classification of B Ds. However, optical/near- infrared 
(IR) s pectroscopy (see, e.g.. iGeballe et al.| [2002: Mc Lean et alj 
120031) has proved impossible for very faint distant BDs with 
current telescopes and very time-consuming for large BD sam- 
ples. To overcome these limitations, new BD classification meth- 
ods have been devised on the basis of near-IR broad-band 
(see, e.g., iLeggett et al.l l20ldl ISchmidt et ail 1201 Ol and refer- 
ences therein) and narrow-band imaging that targets unique 
molecular features of L an d T dwarfs, in particular, water and 
methane bands (see, e. g., iGorlova et al.l 120031 : iMainzer et al.l 
120041: ISpezzi et al.l201 ll) . These techniques are equivalent to ex- 
tremely low-resolution spectroscopy and can be confidently ap- 
plied for statistical purposes, e.g., to detect and classify bona 
fide BDs in large imaging surveys. Our search for T-dwarfs 
in Serpens is based on narrow-band methane (CH4) imaging. 
We target, in particular, the CH4 absorption band centered at 
1.66 /im, which forms in the atmosphere of obje cts cooler than 
Jp.ff $1200 K (defyning the L/T dwarf boundarv: lBurrows et alJ 
2006). This technique has been successfull y applied to a sear ch 
for T-dwarfs in the young IC 348 cluster (Burgess et al. 2009). 

This paper is organized as follows. Section [2] describes the 
WIRCam observations in the Serpens Core and the catalog ex- 
traction procedure. In Sects. [3]|4] we use this catalog to identify 
possible T-dwarfs on the basis of methane absorption bands and 
estimate their intrinsic properties (T e ff> Ay, etc.). In Sect. [3] we 
investigate the nature of our T-dwarfs candidates using comple- 
mentary near to mid-IR data and T-dwarf spectral models. The 
conclusions of this work are given in Sect. [6] 



Table 1. Journal of the observations. 



Field 


Date 

(d/m/y) 


Filter 


T 

(sec) 


Seeing 

(") 


Serpens Core 


17-19/05/2008 


CH 4 off 


2520 


0.6 




17-19/05/2008 


CH 4 on 


7560 


0.6 




20/04/2008 


J 


8405 


0.7 




20/04/2008 


H 


1405 


0.7 




20/04/2008 


K s 


1405 


0.7 



1 This choice is appropriate because we use the distance to Serpens to 
estimate the lower limit to the interstellar extinction towards the cloud 



2. Observations and data reduction 

The imaging observations presented in this work were obtained 
as part of a large CFHT (Canada-France-Hawaii Telescope) key 
program (08AF98) aimed at the characterisation of the low- 
mass population of several young star-forming regions (P.I. J. 
Bouvier). The data were obtained in queue-scheduled observing 
mode between 20 April and 19 May 2008 using the Wide-field 
InfraRed Camera (WIRCam) at CFHT. The WIRCam consists 
of four Hawaii-II2-RG 2048x2048 array det ectors with a pixe l 
scale of 0.3" or 0.15" with microdithering dPuget et alJ l2004). 
The four detectors are arranged in a 2x2 pattern, with a total 
field of view of 21.5'x21.5' and small gaps of 45" between ad- 
jacent chips. The average CCD read-out noise and gain are 30 e- 
and 3.8 e-/ADU, respectively. 

In this paper, we focus on the observations of the Serpens 
core cluster obtained through the CH4off (/lc=l-58 jum, 
AA-O.lfim) and CH 4 on (/t c =1.69 fim, A/}=0.1/mi) filters, along 
with observations through the JHK$ broad-band filters obtained 
using micro-dithering. A single WIRCam pointing centered at 
RA=18 /, 29 m 57i and Dec=+01 fl, 13'"08 5 was sufficient to cover 
the Serpens core cluster. Each tile in each filter was observed us- 
ing a 7-point dithering pattern selected to fill the gaps between 
detectors and accurately subtract the sky background. Exposure 
times, seeing and air mass conditions are summarized in TableQ] 

Individual images were primarly processed using the stan- 
dard CFHT Data Processing and Calibration Pipeline for 
WIRCam (I'iwi; Albert et al., in preparation), which includes 
bias subtraction, flat-fielding, non-linearity correction, cross-talk 
removal, sky subtraction, and astrometric calibration. The astro- 
metric World Coordinate System (WCS) solution is purely linear 
(i.e., assumes a constant scale) and is determined in two steps. 
Knowing the geometry of the detectors (WIRCam has a field dis- 
tortion of about 1 % in the corner of the mosaic), a first fit using 
all 2MASS stars found on the whole mosaic is performed. Then, 
provided enough stars are found, the same fit method operates 
on each individual detectors to refine the WCS solution individ- 
ually. This generally yields an astrometric accuracy of ~1 arc 
second and the RMS scatter of the resulting WCS solution is 
generally about 0.5 arc seconcQ. 

Afterwards, the combination of individual exposures into fi- 
nal stacked images and their photometric calibration were per- 
formed using the SCAMP and SWarp standard packages by 
Terapix, the data red uction centre a t the Institut d' Astrophysique 
de Paris, France (see lMarmoj |2007, and references therein). The 
photometric calibration of the WIRCam data in the broad-band 
filters is done as part of the nominal pipeline reduction using 
2MASS stars in the observed frames and images are renormal- 
ized to an arbitrary photometric zero-point of 30 mag. Typical 
estimated errors in the WIRCam zero point determination are 



(see also Sect.[4]l. 



2 For further details, we defer the reader to the WIRCam webpage: 
http://www. cfht. hawaii. edu/liistrumeiits/Imaging/WIRCam/Iiwi Version! Doc. htmli 
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of the order of 0.05 mag, depending on the number of 2MASS 
stars available, on the variability of those stars and on small 
problems in flat-fielding and/or sky subtraction from night to 
night. A further check/refinement of the extracted JHKs magni- 
tudes was performed using 2MASS photometry for the bright- 
est stars in the field, as is described in the next section (i.e., 
Sect. 12.11 ). Images in the methane filters have no external pho- 
tometric calibration and the CHjoff and CH4011 magnitudes are 
given here on an arbitrary albeit internally consistent scale, so 
that H-CH4on=H-CH4off=0 for unreddened field dwarfs. 

2.1. Photometric catalog 

The source extraction and photo metry were performed by us- 
ing a c ombination o f SExtractor (Ber tin & Arnoutslll996l) and 
PSFex dBertinll201 lb . SExtractor extracts well defined stellar- 
like objects, which are used by PSFex to compute a PSF model 
that is allowed to vary with position on the detector. Then, 
SExtractor uses this PSF model to accurately extract and mea- 
sure the photometry of all the sources detected on the image. 

T-dwarf candidates are expected to be fainter and, hence, 
possibly missed in the CH4011 image because of their CH4 ab- 
sorption. Thus, we first perform the source extraction on the 
CHjoff image and then use this catalog as input list when 
running Sextractor on the CHjon image. To ensure the de- 
tection of all the faint sources present in the images the ex- 
traction criteria used are not too stringent. An object is ex- 
tracted if it complies with the requirement to have three 
contiguous pixels (DETECT JV1INAREA=3) with fluxes 1.5 
times above the background variation (DETECT_THRESH 
=ANALYSIS_THRESH=1.5). In order to minimize the number 
of spurious detections, the SExtractor cleaning efficiency pa- 
rameter was set to the maximum value (CLEAN_PARAM=1). 
An inspection of the images and detections showed that all 
the objects seen by eye are detected. However, this visual in- 
spection showed that a significant fraction of blended sources 
was missed. A higher selection threshold (DETECT_THRESH 
=ANALYSIS_THRESH=5) allowed us to recover those sources 
and we created a final master catalog combining the results of 
the two runs and containing ~130 000 objects. The summary 
of saturation limits and limiting magnitudes, after cleaning for 
saturated objects and obvious spurious detections/artifacts (i.e. 
Sextractor FLAGS <2), is given in Table [2] The approximate 
completeness limit in each filter was derived as the point where 
the histogram of the magnitudes not corrected for extinction 
(Figure [TJ diverge from the dotted line, which represents the 
linear fit to the logarithmic number of objects per magnitude 
bin, calculated over the intervals of better ph otometric accuracy 
dSantiago et al.lll996tlWainscoat et al.lll992l) . 

About 60% of the objects detected in CFLoff and CH4on 
images have a counterpart in the broad-band images. An 
analogous procedure was applied to obtain JHKs pho- 
tometry. We perform the source extraction on the /, H 
and Ks images independently using Sextractor and adopt- 
ing the same configuration parameters as for the methane 
filters, i.e. DETECT_MINAREA= 3 , DETECT.THRESH 
=ANALYSIS_THRESH=1.5 and CLE AN_PARAM= 1 . 
A second run was also performed on the JHKs images 
using an higher selection threshold (DETECT.THRESH 
=ANALYSIS_THRESH=5) to recover the blended sources. 
The final master catalog, obtained by combining the results of 
the two runs, contains ~94000 sources for the / filter, ~81000 
for the H filter and -87000 for the K s filter. Approximately 
4000 of the brightest stars in our field have a counterpart 



Table 2. Saturation limit, limiting magnitudes at the 10cr and 5<x 
level and completeness limit in each filter. 



Filter 


Saturation 


Mag 10tr 


Mag 5ir 


Completeness 




limit 






limit 


CH 4 off 


13.0 


20.1 


20.9 


19.5 


CH 4 on 


13.5 


20.6 


21.3 


19.5 


J 


13.5 


20.8 


21.5 


20.5 


H 


12.5 


19.0 


19.8 


19.0 


Ks 


12.5 


18.6 


19.3 


18.5 



in the 2MASS catalogue. We used this sample to check the 
photometric accuracy and pipeline reduction of the WIRCam 
images and to estimate the difference between the CFHT and 
2MASS photometric systems. The mean magnitude differences 
between the two systems is found to be 0.07, 0.06, and 0.1 
mag for J, H, and Ks, respectively, which is of the order of 
expected zero point uncertainty (~0.05, Sect. [2]). The dispersion 
of the differences is ~0.08 mag form the three filters and reflects 
the possible sources of error mentioned above (i.e., number 
and variability of the available 2MASS stars, non-optimal 
flat-fielding and/or sky subtraction from night to night, etc.; 
Sect. |2). Given that the dispersion is as large as the typical 
magnitude difference between the two systems, we choose to 
not correct our photometry for these offsets. Moreover, this 
comparison hold for the brightest stars in our images, because 
fainter ones are undetected in the 2MASS survey. This implies 
that color effects arising from CFHT filters can not be taken 
into account in this calibration, although they are expected to 
be significant because the WIRCam and 2MASS filters design 
differ substantially. Thus, throughout this work, all the WIRCam 
J, H, and Ks photometry is given in the CFHT Vega system. 
The saturation limits and limiting magnitudes in the broad-band 
filters, after cleaning for saturated objects and obvious spurious 
detections/artifacts (i.e., Sextractor FLAGS<2), are given in 
Tabled The approximate completeness limit in each broad-band 
filter was derived as for the methane filters (Figure[TJ. 

The catalogues from all filters are combined into a single 
database by requiring a positional match better than 1", corre- 
sponding to the astrometric accuracy of our images (see Sect. [2). 



3. T-dwarf candidate selection 



As demonstrated by Burgess et al. (2009), the flux ratio be- 
tween the CH4on and CH4off filters gives a good measure of 
the strength of the methane absorption in T-dwarfs. Specifically, 
the CH4on filter probes the CH4 absorption band centered at 
~1.66 fim and typical of T-dwarfs, while the CHjoff filter covers 
a wavelength range relatively featureless in T-type objects and, 
hence, can be used as the primary conti nuum filter. 

As shown by iBurgess et al.l Q2009) in their Figure 2, the 
CH4on-CH4off color in the CFHT methane filters is roughly 
equal to zero for L-type and earlier type dwarfs and smoothly in- 
creases towards later spectral type, so that T-dwarfs have colors 
above ~0. 1 . Figure |2] shows the CH4on-CH4off against CH4off 
diagram for point-like objects detected in the observed area in 
Serpens core. We use this diagram to select T-dwarf candidates 
following the steps described below: 

- Most of the stars in our field have spectral type earlier than 
L, as demonstrate by the over-density of points at CH4on- 
CTLtoffeO in Figure |2] We use the median CH4on-CH4off 
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14 15 16 17 18 19 12 14 16 18 20 



Fig. 1. Left panels: Photometric errors as a function of PSF magnitudes and relative exponential fit (continuous line) for all the 
point-like sources detected in the CH/toff, CH4on and JHK$ images. Right panels: Number of detection as a function of magnitude. 
The line shows the linear fit to the histogram points used to find the turning point of the distribution, indicating our completeness 
limit. 



color of these stars as a function of the CITtoff magnitude 
to define the reference template with respect to which the 
CH4 absorption feature is sought. The continuous line and 
the line-filled area in Figure|2]show the median locus of field 
stars and their dispersion (cr), respectively. 
- We consider the photometric uncertainty on the ClTton- 
CH 4 off color (6 = y/ACH 4 on 2 + ACH 4 off 2 ) and select a 
first sample of T-dwarf candidates by considering all those 
objects with (CFLton-CFLtofft-tJ color exceeding that of the 
median locus of field stars by more than 3<x. The objects just 



below this limit could still be early T-dwarfs but it would be 
more difficult to extract them because of photometric errors. 
With this criterion ~200 sources were selected. 
- This sample was visually inspected on both the CFUoff and 
CH4on images and the JHKs images. Most of the objects 
turned out to be nebulous detections, ghosts or detector 
cross-talk, spurious detections/artifacts in the vicinity of sat- 
urated objects or image edges. Thus, the initial sample was 
reduced to a shortlist of only 4 point-like candidates. 
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Fig. 2. CH/ion-CHUoff vs. CH401F diagram for the point-like ob- 
jects detected in Serpens core (small grey dots). The continu- 
ous line and the line-filled area show the median locus of field 
stars and their dispersion, respectively. The squares and trian- 
gles show the location of confirmed field L and T dwarfs, re- 
spectively; their magnitudes have been shifted to the distance 
modulus of Serpens. The circles show the position of our the T- 
dwarf candidates; their ID number (see Table [301 is indicated 
and errors bars are within the circle size. The extinction vector 
is shown for Ay=10 ma g. The ClTton-CH/i off vs. spectral type 
calibration relation by Burgess et al. (2009) is also reported on 
the y-axis. 



The PSF photometry of these four T-dwarf candidates in the 
CJ-Uon/CFUoff filters and in JHKs is listed in Table [3] Thumb- 
sized images of the four candidates are shown in Figure[3] 

4. Reddening and spectral type estimates 

Knowledge of the candidates' extinction (Ay) is required in or- 
der to be able to estimate their spectral type and collect further 
information on their nature. 

We do not know a priori whether our candidates are field 
T-dwarfs or young members of Serpens. In both cases, we can- 
not give accurate Ay value for each object based only on pho- 
tometry, because we do not know their distance; they could be 
foreground objects closer than Serpens, cloud members located 
between 26 and 415 pc, i.e . the distance range of Serpens (see 
Table 2 bv iDzib et al.ll2010l) or even being located beyond the 
cloud. Thus, we estimate the range of possible extinction values 
for each candidate as follows: 

- According to the study by Chapman et al. (2009), the mid- 
infrared extinction law in S erpens is consistent with the 
Weingartn er & Draind (1200 ll) Ry= 3.1 diffuse interstellar- 
medium dust model below <1 (Ay <9). We use this 
prescription to define the direction of the extinction vector. 

- Assuming that our candidates belong to Serpens, extinc- 
tion can be estimated using color-color diagrams of CFLjon- 
CtLtoff versus J-H, J-K$ and H-K$ , as plotted in Figure[4] 
The extinction is computed for each candidate using the ex- 
tinction vector and regr essing the objects ba ck towards the 
3 Myr PMS isochrone dChabrier et alj|2000h i.e. the fypi- 
cal age of the Serpens cloud core population dOliveira et alJ 



120091: IWinston etalJl2009h . The extinction values derived 
through this method (Ay My ) are reported in Table[4] 

- R. Gutermuth (private communication) calculated an extinc- 
tion map of the Se rpens core using 2MAS S data and the 
method explained in iGutermuth et all J2007). This map has 
a spatial resolution of 25"x25" and entirely covers the area 
observed by us. By comparing this map with our WIRCam 
images, we estimate that the crowding of our field is such 
that typically ~50 stars fall in the map resolution box. This 
number is suitable for a statistically signif icant estimate of 
Ay in the given box/direction based on the IGuter muth et al. 
(2007) method. Because the extinction map takes into ac- 
count all the extinction occurring along the line of sight, it 
can be used to estimate the maximum extinction towards our 
candidates under the assumption that they are located beyond 
the Serpens cloud. Moreover, the extinction assigned to each 
object (A'™ 1 ' in Table |4| is the average value in a 25" box 
around the object itself and, as such, can be inaccurate for 
very young and highly embedded objects, which present pe- 
culiar extinction due to the material in their envelope/disk. In 
Sect. I5.2l we will see that this is not likely to be the case for 
our candidates. However, the average Ay value computed in 
a given 25" box can be still inaccurate for an object located 
in a specific position within the box, because of the effect 
of small scale differential reddening, presence of cloud fila- 
ments, etc. 

- Finally, to estimate the lower limit to the extinction of our 
candidates, we assumed that they are foreground T-dwarfs 
located within 260 pc of our Sun, i.e. the minimum dis- 
tance estimated for Serpens. The typical diffuse absorption 
is taken t o be 0.7-1 mag/kpc in the Solar Neighborhood 
(see, e.g., Ilchikawa et al.lll982l) . Thus, interloping field ob- 
jects between our Sun and Serpens would be subject to a 
0<Ay <0.2 mag extinction. 

For each candidate, we considered a minimum extinction of 
0<Ay <0.2 mag, the extinction occurring along the line of sight, 
given by the extinction map, and the extinction expected if the 
object belongs to the young population of Serpens (A™ flp and 

Ay M> in Table |4] respectively) and calculated the relative dered- 
dened CH/jon-CH/toff colors. We then used the d-L on-CFUoff 
vs. spectral type calibration bv iBurgess et a l. (2009) (see their 
Figure 2) to estimate the spectral type range consistent with 
these dereddened co lors. Note that the calibration relation by 
Burg ess et alJ (120091) is based on field dwarfs; the actual col- 
ors of our LT dwarfs might be slightly different because of 
the effects that redu ced gravity has on the opacities (see, e.g., 
ICasewell et al.ll20lll and references therein). 

The spectral type estimates for our candidates are reported 
in Table|4] together with the corresponding ef fective temperatur e 
range according to the calibration relation bv lVrba et al.l J2004). 

5. On the nature of the T-dwarf candidates: are they 
young members of Serpens? 

In this section we use color-magnitude diagrams (CMDs), color- 
color (CC) diagrams and spectral energy distributions (SEDs) 
to assess whether the optical/IR properties of our T-dwarf can- 
didates are consistent with those of field objects or indicate the 
presence of IR excess emission typical of very young objects, 
which would in turn indicate membership to the Serpens core 
cluster. This analysis is performed using our WIRCam data in 
combination with r and /-band imaging data from MegaCam at 
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Table 3. Photometry of the T-dwarf candidates. CH4011 and CH4off magnitudes are in the WIRCam system, while riJHKs photom- 
etry is given in the Vega system. 



ID 


R.A. J2000 


Dec. J2000 


r CFHT 


Is ubaru 




CH 4 off 


CH 4 on 


J 


H 




K s 




(hh:mm:ss) 


(dd:mm:ss) 




















1 


18:29:56.58 


+01:18:35.93 




22.98d 


=0.03 


17.77+0.05 


18.07+0.05 


19.10+0.05 


17.71 = 


t0.06 


16.99+0.05 


2 


18:29:57.64 


+01:19:40.61 




22.67d 


=0.03 


17.61+0.05 


18.10+0.05 


18.84+0.05 


17.64= 


t0.06 


17.07+0.06 


3 


18:30:27.89 


+01:14:52.22 




23.98=i 


=0.10 


18.75+0.05 


19.16+0.06 


19.38+0.06 


19.08= 


t0.16 


18.90+0.27 


4 


18:30:37.24 


+01:18:37.68 


23.57+0.07 


22.14+0.02 


18.61+0.06 


18.84+0.06 


19.71+0.07 


18.67= 


t0.09 


18.25+0.10 




Fig. 3. Images of the four T-dwarfs candidates in each of the WIRCam@CFHT filter used in this work, as indicated in the labels. 
Each snapshot covers an area of 18"x 15"; north is up and east to the left. 



Table 4. Estimates of visual extinction, spectral type and effec- 
tive temperature for our T-dwarf candidates based on their dered- 
dened CH4on-CH4off colors. Ay My is the visual extinction when 
assuming that the object belongs to the 3 My old population of 
Serpens, while A™" 1 ' is derived from the extinction map by R. 
Gutermuth (private communication). 



ID 


A map 
A V 


A iMy 


Spec. Type 


T f // 




(mag) 


(mag) 




(K) 


1 


9.2 


9.4 


T2-T4 


1390-1290 


2 


7.9 


7.1 


T4-T5 


1290-1190 


3 


4.2 


1.0 


T3-T4 


1360-1290 


4 


2.8 


5.5 


T1-T3 


1390-1360 



CFHT and Suprime-Cam at the Subaru telescope, and Spitzer 
imaging. 



5. 1 . WIRCam near-IR color-magnitude and color-color 
diagrams 

Figure [5] shows the position of our candidates on the J vs. J - 
K s CMD and the J - H vs. H - K s CC diagram before and 
after correction for interstellar reddening. In order to draw the 
de-reddening vector for our candidates, we adopt the maximum 
extinction expected towards our candidates (i.e., the maximum 
between the A'y' p and A^, Ml value from Table|4|; in this way, we 
display the maximum possible shift of their position in the CMD 
and CC diagram due to reddening effects. 

We also plotted the 5 Gyr DUSTY an d COND models 
and the confirmed field L and T dwarf s (IChiu et al .1 120061: 
iGolimowski et alj|2004t iKnapp et alj|2004l) to highlight the po- 
sition expected for the field dwarf sequence. The locus ex- 
pected for the young members of Serpens is indicated by 
the 3 Myr DUSTY and COND models and the young stel- 
lar object (YSO) population confirmed by the Spitzer core 
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Fig. 4. CH4on-CH4off vs. J - H, J - K and H - K for the point-like objects detected in Serpens core (small grey dots). The circles 
with error bars are the T-dwarf candidates in thi s study; their ID number is indicated. The squares and triang les show the location 
of confirmed field L and T dwarfs, respectively (Chiu et al. 2006; Golimowski et al. 2004; Knapp et al. 2004). The position of the 
CFHTJ0344+3206 mid-T dwarf candidate (green diamond) is plotted for comparison. The lines show the 5 Gyr and 3 Myr COND 
isochrones (dashed lines) and the 5 Gyr and 3 My DUSTY isochrones (continuous line), as indicated in the labels. The extinction 
vector is shown for Ay=10 mag. 



to disk (c2d) survey dHarvev et al.ll2007t I Winston et aT]|2007l) . 
Finally, we plot f or comparison the fa i nt mi d-T type object 
S Ori 70 found by Zapa tero Osorio et al.l (120021) towards the di- 
rection of the young cr Orionis cluster and the dereddened posi- 
tion of CFHT J0344 +3206, a mid-T dwarf candidate found by 
iBurgess et al] j2009h in the direction of IC 348. The DUSTY 
and COND models and the positions of the field LT-dwarfs, 
S Ori 70 and CFHTJ0344+3206 are shifted to the minimum 
distance estimated for Serpens (i.e. 260 pc). Note that it is still 
unclear whether S Ori 70 is a field brown dwarf or a young plane- 
tary mass member of cr Orionis; however, proper motion and the 
near- and mid-IR colors measurements by Zapatero O sorio et alj 
(2008) support its membership to the cluster, with an estimated 
mass in the interval 2-7 M . 

The position in both diagrams of two of our candidates, 
namely ID 1 and 2, is consistent with them being mid-T dwarfs. 
However, their dereddened magnitudes are brighter than ex- 
pected for young T-dwarfs at the distance of Serpens and, hence, 
they might be field objects. They suffer from significant inter- 
stellar extinction (Table [U an d therefore can n ot be foreground 
objects. As discussed bv TDelorme eTail (|2008), the only back- 
ground contaminants may eventually be high redshift quasi- 
stellar objects (QSO; z>6) or heavily reddened low redshift 
QSOs, which appear point-like and share the same very red col- 
ors as cool dwarfs. Moreover, high redshift/heavily veiled QSOs 
can present a moderate methane-like absorption signal in the 
CHton-CHjoff difference du e to the combination s of specific 
redshifts and emission lines (Richards et al. 2003). In Sect. 15.21 
we use additional imaging in the i and Spitzer pass-bands to as- 
sess whether this is the case for ID 1 and 2. 

ID 3 does not present a significant interstellar extinction 
and its near-IR colors are consistent with it being a field mid-T 
dwarf. Its position on the J vs. J - K$ diagram and spectral type 
are very similar to those of the mid-T type dwarf S Ori 70 in the 
cr Orionis cluster and the T-dwarf candidate CFHTJ0344+3206 
in the direction of IC 348. However, it is not redder than the se- 
quence defined by field T-type dwarfs in the /-//vs. H-K$ CC 
diagram, as young objects should be. Indeed, the K-band flux 
is sensitive to gravity and it is fainter for young objects, espe- 
cially for late T-dwarfs, because of the strongly reduced gravity 



( Zapat ero Osorio et al. I I20081) . Since it is brighter than field ob- 
jects at the distance of Serpens (Figure [5] left panel), ID 3 may 
be a foreground field T-dwarf. 

The interpretation of the CMD and CC diagram for candi- 
date ID 4 is controversial, because it lies in both diagrams in a 
"transition region" mainly populated by field dwarfs but where 
faint and/or highly veiled YSOs might still be found. It cannot be 
a foreground object as it is significantly extincted (Tabled and 
has near-IR dereddened colors consistent with it being an early 
T-dwarf. It is brighter than field objects (Figure [5] left panel), 
but this is still consistent with it being young (the larger radius 
of YSOs make them brighter than older objects with the same 
spectral type). Thus, ID 4 may be a cluster member. 

It is worth noting that these very same conclusions can be 
drawn from the inspection of the CC diagrams involving colors 
in the CH4 filters, which we presented in Figure[4] 



5.2. Complementary MegaCam, Suprime-Cam and Spitzer 
data 

The Serpens core cluster wa s observed on 2010 -06-11 in the r 
and i bands with MegaCam (Boulade et al. 2003 ) at CFHT and 
on 20 08-06-05 in the /-band with Suprime-Cam (Mivazaki et al.l 
2002), the wide field camera located at the prime focus of the 
Subaru Telescope. For details about these observations and the 
relative data reduction and calibration, we defer the reader to 
Bouy at al. (2012, in preparation). In Table [3] we report the ri 
photometry of our four T-dwarf candidates calibrated to the Vega 
system. Only one of our candidates (ID 4) was detected in both 
filters, while the remaining three objects (ID 1, 2 and 3) were 
detected only in the /-band. 

We also searched for possible mid- to far-I R flux measure- 
ments of our T- dwarf candidates in the AKARI ( Yamauch i et al.l 
120111) . WISE dWright et all 1201 Ol) and Spitzer databases. The 
AKARI and WISE all-sky surveys turn out to be too shallow 
for the detection of our T-dwarf candidates. However, three of 
them (ID 1, ID 2 and ID 3) have been detected in the deeper 
Spitzer c2d Survey, which mapped Serpens together with other 
five nearby molecular clouds, usi ng both IRAC and MIPS, the 
two imaging cameras on Spitzer dEvans et al.ll2009|) . ID 4, the 
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2 4 0.0 0.5 1.0 1.5 

J - K s H - K s 



Fig. 5. J vs. J - Ks (left panel) and J - H vs. H - Ks diagrams (right panel) for the point-like objects detected in Serpens core 
(small grey dots). The circles with error bars are the T-dwarf candidates in this study; their ID number is indicated and the arrows 
indicate their position after correction for interstellar reddening. The Serpens core class II (black dots) and class III (asterisks) YSO 
population, the S Ori 70 mid-T dwarf (sky-blue diamond) and the CFHTJ0344+3206 mid-T dwarf candidate (green diamond) 
are plotted for comparison. The squares a nd triangles show the location of confirmed field L and T dwarfs dChiu et al. 2006; 
Golimow ski et al.l l2004; Kna pp et"al] |2"004). respectively, shifted to the distance of Serpens. Lines show the 5 Gyr and 3 My COND 
isochrones (dashed lines) and the 5 Gyr and 3 My DUSTY isochrones (continuous line). 



faintest of our candidates, was not detected by the Spitzer c2d 
Survey. The IRAC and MIPS 24//m fluxes for ID 1, ID 2 and 
ID 3 are reported in Table [5] We also computed for these three 
candidates the SED slope between the K-band (2.2/mi) and the 
MIPS band at 24pm (a[K&24^m], Table ID, which is traditio nally 
used to assess the IR class of YSOs (Lada & Wilkingll 984). The 
a[K&24iim] we obtain are upper limits to the actual SED slopes, 
because only flux upper limits are available at 24pm for our ob- 
jects. 

In Figure [6] (left panel) we compare the position of our 
candidates in the Ks - [3.6] vs. [3 .6]-[4.5] CC diag ram with 
the locus expected for field dwarfs dPatten et al. 2006), t he se- 
quence defined by YSOs belonging to Serpens dHarvev et alj 
l2007t IWinston et all 120071) and the position of S Ori 70 and 
CFHTJ0344+3206. 

In Figure [6] (right panel) we use the i - J vs. J-[3.6] CC 
diagram to compare the colors of our candidates with the lo- 
cus of YSOs in Serpens, the position of S Ori 70 and the colors 
expected for background QS Os, which might contaminate our 
sample. As demonstrated by Bou v et aTl (12009). very low mass 
YSOs (with or without mid-IR excess associated to the presence 
of a circumstellar disc) and QSOs occupy two very distinct areas 
of this diagram, thus it can be used to effectively separated the 
two populations. 

Finally, in Figure [7] we present the SEDs of our four can- 
didates from the optical up to 24/im. The line-filled area in 
each SED plot represents the flux variation depending on the 
adopted reddening correction (see Sect. [4] and Table @J. The 
sp ectrum overplotted to e ach SED is the CON D model spectrum 
bv lAllard et all d2001l) and lBaraffe etall d2003l) with the same ef- 
fective temperature as the object, normalized to the J-band flux; 



this model represents the pure stellar photospheric emission ex- 
pected for an object of the given T e yy. 

The Spitzer colors of ID 1 and ID 2 (Figure [6) are in gen- 
eral agreement with those of class III objects with thin or no 
disks. This is also in agreement with the shape of their SEDs, 
which might present a small IR excess emission only at 24 pm 
(Figure |7]i. For both objects we find a[K&24^m] ^ -1, which 
is typical of class III objects with no pr ominent IR excess 
dGreene et alj Il994t lLada & Wilkind fl984). Thus, the overall 
Spitzer colors/SED for these two objects indicate that they do not 
possess a clear IR excess emission typical of very young objects. 
Their position in the i - J vs. 7-[3.6] diagram indicates that their 
colors, when corrected for the estimated interstellar extinction 
(Tabled, approach the locus expected for QSOs. Since they are 
too bright to be high-redshift QSOs, the most likely possibility 
is that they are heavily reddened low redshift QSOs. Indeed, we 
estimated that the Ha broad emission line of QSOs could con- 
taminate the magnit ude in the methane filters at redfshift z~ 1 .4 
dRichards et alJl2003h . 

The SED of ID 3 (Figure |7]i might present IR excess start- 
ing from 8yum with a slope (a[K&2Aim\ =0.27+0.04) typical of a 
class II YSO with a thick disk. However, only flux upper lim- 
its are available for A >5.8/mi and, hence, this (X[K&2Afim] value 
is not robust. The interpretation of the Ks - [3.6] vs. [3.6]-[4.5] 
CC diagram (Figure |6]l is very uncertain because of the big er- 
ror bars. However, the position of ID 3 appears to be consistent 
with the locus of field T-dwarfs. In particular, its [3.6] - [4.5] 
color is bluer than expected for young mid-T dwarfs such as 
S Ori 70 and CFHTJ0344+3206. However, assuming that a 
thick disk surrounds this object, its [3.6] - [4.5] color, peculiar 
for a young objects, might still be explained in term of different 
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disk inclination with respect to S Ori 70 and CFH TJ0344+3206 
Indeed , as shown in the top panels of Figure 7 by Robitaille et alJ 
(2006), disk inclination significantly alter the SED shape at mid- 
IR wavelengths (3-10 pm) for very low-mass objects (<0.2 M ). 
Finally, its position in the i - J vs. 7-[3.6] diagram is clearly 
consistent with the YSO locus and, in particular, it is very close 
to the location of S Ori 70. More robust imaging data at mid-IR 
wavelengths and/or follow-up spectroscopy is needed to draw a 
firm conclusion about this object. 

ID 4 is not detected in any of the Spitzer pass-bands 
and, hence, we can not assess the presence of a circumstel- 
lar disk around this object. For the sake of completeness, the 
optical/near-IR SED of ID 4 is shown in Figure|7] The observed 
i — J color of ID 4 (i.e., 2.43) i s bluer than expected for early/mid 
T-dwarfs (see, e.g., Table 3 bv lHawlev et alJl2002h . while its ob- 
served and dereddened r - i color (1.43 and 0.82, respectively, 
adopting the Serpens extinction map (A'y' p in Table |4]l are con- 
sistent with the r - i colors expected for early/mid T-dwarfs (see, 
e.g., Table 3 bv lHawlev et al.ll2002l) . in agreement with its near- 
IR colors (Figure |5) . Thus, our conclusions about this object re- 
main as in Sect. 15.11 ID 4 may be a cluster member, although 
our data are not sufficient to conclusively prove this and spectro- 
scopic follow-up is required. 



6. Summary and conclusions 

We presented a deep methane imaging survey for planetary-mass 
T-dwarfs in the Serpens Core cluster. Our survey covers a field 
of view of 21.5'x21.5' and about 127000 sources were detected 
in both the CFLjon and CF^off filters. 

We identified four potential T-dwarfs from their 1.6jum 
methane absorption band and used complementary riJHKg 
broad-band imaging and mid-IR flux measurements from the 
Spitzer c2d Survey to investigate their stellar and disk proper- 
ties. 

Two of our candidates (ID 1 and ID 2) suffer from significant 
interstellar extinction and have near-IR colors similar to mid-T 
dwarfs. However, they are too bright to be planetary-mass mem- 
bers of Serpens and do not present the typical IR excess emission 
expected for young objects surrounded by a disk. Our analysis 
indicates that their dereddened near-IR colors approach the val- 
ues expected for QSOs. In particular, since they are too bright 
to be high-redshift QSOs, they could be heavily reddened low 
redshift QSOs. 

Our candidate ID 3 does not present a significant interstel- 
lar extinction and its near-IR colors are more consistent with it 
being a field mid-T dwarf. It is brighter than field objects at the 
distance of Serpens and, hence, it may be a foreground T-dwarf. 
It is detected by the Spitzer c2d survey but only flux upper lim- 
its are available for A >5.8pm and, hence, we can not assess 
the presence of a possible disk around this object. However, it is 
worth to keep this possibility open, because the position on the 
J vs. J - Ks and i - J vs. 7-[3.6] diagrams and the estimated 
spectral type of ID 3 are very similar to those of other young 
T-dwarf candidates, i.e. S Ori 70 and CFHTJ0344+3206 in the 
direction of IC 348. 

Finally, ID 4 is significantly extincted, which exclude the 
possibility of it being a foreground object, and has dered- 
dened optical and near-IR colors mainly consistent with those of 
early/mid T-dwarfs. It is brighter than field objects and, hence, it 
might be very young, as YSOs have larger radii and are brighter 
than older field objects of the same spectral class. Thus, ID 4 is 
a promising young T-dwarf candidate in our sample. 



Considering the magnitudes of our candidates ID 3 and 4 
(Table |3]l and their T e ff (Table |4]l based on the dereddened 
CH4on-CH4off colors, we estimated their masses assuming that 
they belong to the Serpens core cluster (i.e., age~3 Myrs and 
distance in the range 260-415 pc) a nd using both the C OND and 
DUSTY evolutionary models bv lChabrier et alJ d2000t) . The es- 
timated mass is in the range 2-4 Jupiter masses for both ID 3 
and 4. Therefore, if they truly belong to the Serpens Core clus- 
ter, they would be amongst the youngest, lowest mass objects 
detected in a star-forming region so far. 

Follow-up spectroscopy is required to confirm the spectral 
type and reddening of our four T-dwarf candidates and, hence, 
draw firm conclusions about their nature. Additional deep imag- 
ing at mid-IR wavelengths is needed to clarify the presence of a 
possible disk around ID 3. 
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Fig. 6. Left panel: K$ - [3.6] vs. [3.6]-[4.5] color-color diagram. The circles with error bars are our the T-dwarf candidates; their ID 
number (see Table[3]|4| is indicated and the arrows indicate their position after correction for interstellar reddening. The Serpens core 
class II (black dots) and class III (asterisks) YSO population, the S Ori 70 mid- T dwarf (sky - blue d iamond), the CFHT J0344+3206 
mid-T dwarf candidate (green diamond) and the field dwarf sequence from iPatten et"a"D (120061) (line-filled area) are plotted for 
comparison. Righ panel: i - J vs. 7-[3.6] color -color diagram. Symbols ar e as in the left panel. The dot-dashed ellipse marks the 
area occupied by QSOs in the SWIRE catalog (Hatziminaoglou et al. 2008). 



Table 5. Spitzer fluxes and SED slopes («[x&24^m]) f° r the three T-dwarf candidates in this study detected by the Spitzer c2d survey. 



ID 


c2d Name 




IRAC 3.6pm 


IRAC 4.5pm 


IRAC 5.8pm 


IRAC 8pm 


MIPS 24pm 








(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


1 


SSTc2dJ182956.6+011836 


-0.88+0.02 


0.0869±0.0084 


0.0610+0.0065 


0.0930±0.0307 


0.0464+0.0405 


< 0.1410 


2 


SSTc2dJ182957.6+011941 


-1.06+0.02 


0.0667±0.0084 


0.0476±0.0102 


0.0393±0.0403 


<0.0028 


<0.0836 


3 


SSTc2dJ183027.9+011452 


0.27+0.04 


0.0295±0.0078 


0.0121+0.0108 


<0.0060 


<0.0590 


< 0.3930 
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Fig. 7. Spectral energy distributions of the T-dwarf candidates in this study. Flux errors are smaller than the symbol size and 
flux upper limits are marked with an arrow. The line-filled ar ea represents the flu x va riation depending o n the adopted reddening 
correction. The oveplotted spectrum is the COND model by Alla rd et al.l d2001h and iBaraffe et alj d2003l) with the same T e ff as 
estimated for each object (indicated in the label and in Table HJi normalized to the /-band flux. 
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